
HREC-l 156-1 
LMSC/HREC D] 62657 


LOCKHEED MISSILES A SPACE COMPANY 
HUNTSVILLE RESEARCH A ENGINEERING CENTER 
HUNTSVILLE RESEARCH PARK 
4800 BRADFORD DRIVE, HUNTSVILLE, ALABAMA 


SPACE SHUTTLE VEHICLE ROCKET 
PLUME IMPINGEMENT STUDY 
FOR SEPARATION ANALYSIS 


TASKS II and HI 


DEFINITION AND PRELIMINARY PLUME 
IMPINGEMENT ANALYSIS FOR THE MSC BOOSTER 


November 1970 


Contract NAS9-11156 


by 


Carl J. Wojcie chow ski 
Morris M. Penny 
Robert J. Prozan 


APPROVED: 



George D. Reny, Mapraj 


ager 

Aeromechanics Dept. 




PftSCLTlMG r.\ 


p-lank not rn^im 


LMSC/HREC D162657 


FOREWORD 

This report presents the results of Tasks II and III per- 
formed by Lockheed's Huntsville Research & Engineering Center 
for the Flight Technology Branch of NASA's Manned Spacecraft 
Center, Houston, Texas. This work was performed under 
Contract NAS9-11156 and was monitored by Mrs. I. H. Fossler. 

Task I was previously issued as Lockheed dccument 
LMSC/HREC D162574, dated August 1970. 
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SUMMARY 

This document presents the results of a space shuttle plume impingement 
study for the Manned Spacecraft Center configuration. This study was con- 
ducted as two tasks which were to: (1) define the orbiter main stage engine 

exhaust plume flow field, and (2) define the plume impingement heating, force 
and resulting moment environments on the booster during the staging maneu- 
ver. To adequately define these environments during the staging maneuver and 
allow for deviation from the nominal separation trajectory, a multitude of 
relative orbiter/booster positions are analyzed which map the region that con- 
tains the separation trajectories. The data presented in this report can be 
used to determine a separation trajectory which will result in acceptable im- 
pingement heating rates, forces, and the resulting moments. The data, pre- 
sented in graphical form, include the effect of roll, pitch and yaw maneuvers 
for the booster. 

Quasi-steady state analysis methods were used with the orbiter engine 
operating at full xhrust. To obtain partial thrust results, simple ratio 
equations are presented. Using a simplified thrust buildup curve, the heat- 
ing rate lime histories to selected points on the booster vertical tail are used 
in a ur.e-dimensional heat conduction program to calculate the surface 
tempe rature -time histories. These surface temperature -time histories are 
calculated for two separation trajectories. The maximum calculated surface 
temperature at the leading edge of the vertical tail is 1800°F for both of the 
separation trajectories analyzed. The titanium skin on the sides of the verti- 
cal tail reached 480°F and 440°F for the two separation trajectories analyzed. 
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The plume impingement force and moment analysis for the booster dis- 
closes several cases which could be potential problems during staging. The 
analysis revealed that if the booster was perturbed in the yaw plane during 
the staging sequence, the orbiter plume will produce am unstable condition for 
the booster, i.e., moments generated during plume impingement will augment 
the booster yaw motion. The analysis also reveals that the orbiter plume im- 
pingement will produce moments in pitch that will augment any positive pitch 
motion of the booster. 
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NOMENCLATURE 

area, units consistent with surface equations 
engine area ratio, dimensionless 

stagnation pressure coefficient behind a normal shock, dimensionless 
error function of the molecular speed ratio, Sr 

denotes unit vector quantity 

force vector, F i + F U F k, lb- 
x y J z f 

Knudsen number, dimensionless 

direction cosines between the subshape axis system and 
the local force direction 

number of elemental areas used to describe a given subshape 

A a A 

moment vector, M i + M j + M k, ft-lb, 

’ x y* z f 

component of Mach number normal to the elemental area, dimen- 
sionless 

number of subshapes, used in Eqs.(C-6) and (C-8) 

engine propellant oxidizer to fuel ratio, dimensionless 

local static pressure, psi 

body/local impact pressure, psi 
2 / 2/0 

gas constant, ft /sec / R 
molecular speed ratio, v/ 2RT 

distance vector, units consistent with subshape surface equation 
local static temperature, °R 
incident molecular temperature, °R 

reflected molecular temperature and is equal to the body 
temperature for complete energy accommodation 

transformation matrix from a given subshape coordinate system 
to the booster reference coordinate system 
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Greek 
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p, y , t) 
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entropy, ft /sec / R 
velocity, ft/ sec 

plume Cartesian coordinate system 
plume cylindrical coordinate system 


local impingement angle of velocity vector with respect to the 
local body surface tangent (degrees) 

1/2 RT 

ratio of specific heats, dimensionless 

local surface unit normal vector written in local subshape co- 
ordinate system 

density of molecules 

density, slugs/ft 

direction cosines between the subshape local axis system and 
the mass velocity 

indicates summation 

unit vector 


Subscripts 

oo 

tran 

cont 

f.m. 

B.N. 

e 

exit 


indicates free stream conditions 
indicates transitional flow regime 
indicates continuum flow regime 
indicates free -molecular flow regime 
indicates booster nose 
engine lip angle, degrees 
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Section 1 
INTRODUCTION 

The separation of launch vehicle stages during boost usually results 
in the rocket exhaust plume from the upper stage impinging on the lower 
stage. For expendable boosters, the impingement has historically been of 
minor consequence. For the fully reusable space shuttle vehicle, however, 
the impingement problem is significant because control of the lower stage 
must be maintained, and the structure must not be oamaged. The impinge- 
ment environment imposed on the booster can be severe if ignition of the 
orbiter main propulsion system occurs while the stages are close together. 

If ignition of the orbiter is delayed until after the booster is far away, result- 
ing in the impingement environment being either inconsequential or non- 
existent, then payload performance penalties can occur. An assessment is 
needed of the significance of the staging plume impingement problem along 
with the definition of an "optimum" staging trajectory. This optimum tra- 
jectory should provide an acceptable tradeoff between impact of the vehicle 
design to accept the plume impingement environment and penalties incurred 
in system performance and/or operation. 

Task II of this contract defines the engine plume flow field of the Manned 
Spacecraft Center (MSC) Space Shuttle o rbite r and is presented in this docu- 
ment as Appendix A. Task III, presented in the following sections, defines 
the orbiter Dlume impingement heating, force, and moment environments on 
the MSC booster during the staging maneuver. 
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1.1 PROBLEM DESCRIPTION 

The geometrical orientation for the MSC Space Shuttle arrangement is 
shown in Fig. 1.1. A reference coordi- dte system is located at the nose of 
the booster on the vehicle centerline. All forces and moments presented are 
referenced to this point. A simulated multiple nozzle was used to represent 
the main stage engines of the orbiter. The simulated nozzle has the same 
mass flow rate as the orbiter engine cluster and approximately the same plume 
shape. Details of the orbiter plume generation are presented in Appendix A. 

1.2 DEFINITION OF THE PLUM1 'PINGEMENT CHARACTERISTICS 

The impingement characteristics were calculated by Lockheed's Plume 
Impingement Computer Program. Basically, the procedure for obtaining the 
impingement characteristics is as follows: 

The vehicle (composite shape) and the particular motor are located with 
respect to the reference coordinate system. The composite shape is represented 
by a number of geometrically simple subshapes (each of which has its own local 
coordinate system). These subshapes are then subdivided into elemental areas. 
Loc? I flow conditions on a given body are obtained by locating fhe centroid of 
each elemental area in the free plume from a search of the plume flow field 
which is stored on magnetic tape. The flow properties thus obtained are 
assumed to act uniformly over the elemental area. 

The body local flow regimes were determined by the value of the local 
Knudsen number. In this analysis, continuum flow was assumed to exist for 
Knudsen numbers less than 0.01; transitional flow exists for Knudsen numbers 
between 0.01 and 10.0; and Knudsen numbers greater than 10.0 define the 
"free molecular flow" regime. The characteristics lengths used in the 
Knudsen number calculations were a representative dimension (element area 
dimension) of the body in question so that the elemental area local flow regime 
is determined by comparing the mean free molecular path of the local plume 
to the body dimension, 
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Plume 

-Boundary 

Shock Wave 



Drawing not to scale 

Force and moment arrows 
indicate positive directions 

Booster reference coordinate 
system located at booster nose. 


Fig. 1-1 - Three-Dimensional Perspective of the Booster Reference 

Coordinate System and the Orbiter Plume Coordinate System 
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The plume coordinates with respect to the sujface immersed in the 
plume are shown in the sketch below: 



The body local flow properties are computed for a surface immersed in a 
free plume. The properties for a point r, x in the plume coordinate system 
are obtained from a search of the flowfield data stored on magnetic tape. 

The tape search routines use a four-point interpolation scheme to obtain the 
flow properties at the point (r, x). 
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The plume velocity vector is given by: 


where 


A A 

v i + vi, 

x r x 


( 1 ) 


v = v cos e 

X 

V = V sin0 
r 


( 2 ) 


The velocity vector at the body point with respect to the radial orientation 
within the plume is given by: 


where 


A A A 

V = V i + V j + V k 
x y J z 


V =V cost/?, cosio-y/i 
y r 


V = V sine/?, sin</> = z/r 
z r 


(3) 


( 4 ) 


A 

e 

n 



= inclination of the local velocity 
vector with respect to the local 
surface tangent 


The body local impingement angle (shown in the sketch above) is then given 
by 

a = 90.0 - 57.295778 (cos' 1 (e . e)). (5) 

' ' v n ' ' 
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Section 2 

DISCUSSION OF RESULTS 

The MSC booster configuration analyzed is depicted in Fig. 2 -la. The 
overall dimensions used are shown as well as the locations of each of the 
vehicle subshapes (i.e., wings, tail, fuselage, etc.) analyzed. The solid 
outlines show the vehicle shape as analyzed while the dashed outlines show 
the actual vehicle shape. Fig. 2 -lb shows a three-dimensional perspective 
view of the booster reference coordinate system used as well as the orbiter 
plume reference coordinate system used. 

Shown on Fig. 2-lc are the assumed directions of positive yaw, }j/, 
pitch, 0 and roll, <p , of the booster as used to generate the data in Figs. 2-2 
through 2-22. Also indicated is the coordinate transformation matrix used 
to locate the booster nose with respect to the orbiter plume coordinate 
system (Fig. 2-lc). It should be noted that the order of axis rotation to ob- 
tain the transformation matrix is the yaw, pitch and roll directions, respec- 
tively. A discussion of the method used to obtain the properties at a point on 
the booster with respect to the orbiter plume is given in Appendixes B and C, 
respectively. 

In the plume impingement analysis, the various components of the 
booster were represented by four basic subshapes (Fig. 2- Id). These were: 
(1) a cylinder for the upper portion of the fuselage; (2) flat plates for the 
lower sides of the fuselage; (3) a conical nose; and (4) airfoil shape NACA 
0012-64 for the wing, vertical and horizontal tails, respectively. A detailed 
description of the airfoil subroutine is presented in Appendix B. 

In the numerical integration scheme, the cylindrical portion of the 
fuselage and the conical nose section were divided in rings of 100 inches in 
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increment with each ring further subdivided into 10-degree radial increments. 
All wing sections were divided into approximately 50 -inch increments in the 
span direction. As discussed in Appendix B, a cosine distribution function 
is used to distribute the chord direction. This scheme biases the point dis- 
tribution to wing-leading edge where largest gradients in surface curvature 
occur. Twenty increments in the chord direction were used on the wing, 
horizontal and vertical tail surfaces. Uniform flow was assumed to exist 
over each elemental area of the respective surfaces with vacuum conditions 
considered to exist on the surfaces not exposed to the engine exhaust plume. 
Figs. 2-le through 2-lm show various subshapes used to represent the 
booster. All wing surfaces (see Appendix B) were represented by airfoil 
shape NACA 0012-64. 

A graphic representation of the orbiter engine exit plane/booster sepa- 
ration distances used to map the plume impingement environment on the booster 
are shown in Fig. 2 -In as "o" points. A plume overlay is provided with 
Fig. 2 -In to aid in visualization of the orbiter exhaust plume impingement on 
the various booster surfaces. The basic plume used in the analysis (see 
Appendix A) was generated assuming the staging altitude to be 283,000 feet 
with a freestream Mach number of 12.08. The resulting plume shape and dis- 
tribution of plume flow field properties are shown on Figs. A-l through A-6. 

2.1 IMPINGEMENT FORCE AND MOMENT ANALYSIS 

The resultant forces and moments acting on the booster resulting from 
the orbiter main stage engine exhaust plume impingement are presented in 
Figs. 2-2 through 2-22. The data presented comprise a map of the expected 
forces and moments acting on the booster over the region of the separation 
trajectory of the orbiter vehicle (Fig. 2- In). For this study the reference 
coordinate system for the moments is located at the nose of the booster 


2-2 


LOCKHEED • HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC/HREC D162657 


(Fig. 2- lb). The moment data as presented can be referenced to the booster 
center of gravity by the following relation: 



g- 




g- 


x F 


b.n. 


(3.1) 


where the vector 
booster c. g. 


quantity r c 


g- 


is the vector from the booster nose to the 


The impingement data presented were generated assuming the orbiter 
main stage engine to be operating at full thrust. However, the force and 
moment components at any time, t, during the engine buildup to full thrust 
may be approximated by multiplying the force and moment data by the ratio 
of the chamber pressure at the time, t, to an engine chamber pressure of 
3000 psia (see discussion in Appendix C). 

Figs. 2-2 through 2-4 show the impingement force and moment maps for 
the zero pitch, roll and yaw cases. All physical dimensions are with respect 
to the reference coordinate system and are located at the booster nose and indi- 
cate various positions of the orbiter engine exit plane with respect to the 
booster nose. Since that portion of the booster which is subject to orbiter 
exhaust plume impingement is symmetrical about the booster centerline 
(Fig. 2- lb), the moments acting about the booster x (roll) and z (yaw) axes 
are zero with a zero force component in the y direction also. Due to the 
large change in the impingement forces and the resulting moments over the 
region of interest a semi-log interpolation scheme should be employed when 
using the data in tabulated form. The very large increase in the magnitude of 
the impingement data is due to the booster surfaces becoming more fully 
immersed in the orbiter exhaust plume as the separation distance between the 
vehicle increases. 


Figs. 2-5 through 2-7 present impingement data for the booster due to 
a 5-degree pitch maneuver of the booster 'ith respect to the orbiter exhaust 
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plume. This maneuver represents a rotation of the booster nose into the 
orbiter exhaust plume. Figs. 2-8 through 2-10 present the impingement 
data for the booster pitched minus 5-degrees with respect to the orbiter 
exhaust plume. An examination of the figures which show the positive mo- 
ment direction and the impingement data of Figs. 2-7 and 2-8 indicates the 
following interesting problem. When the booster is rotated in a positive 
pitch maneuver (+ pitch) the resulting moment will tend to augment the rota- 
tion (continue to rotate the vehicle in the positive pitch direction). However, 
the resulting pitching moment for a maneuver which tends to pitch the 
booster vehicle in the negative pitch direction (-pitch) will tend to counter- 
act the original direction of rotation; i.e., the direction of rotation due to 
the pitching moment would be opposite the original direction of rotation of 
the booster nose. The above discussion is true for all cases in which the 
plume impingement center of pressure, (booster location where the resultant 
plume impingement force is acting) is aft of the booster center of gravity. 
The reverse is true when the plume impingement center of pressure moves 
forward of the booster center of gravity (which would occur at large separa- 
tion distances). 

Figs. 2-11 through 2-16 present the impingement data for a 5 -degree 
roll maneuver of the booster while Figs. 2-17 through 2-22 present the irr - 
pingement data for a 5 -degree yawing maneuver. Two .ems of consequence 
should be noted here. First, similar to the positive pitch case discussed 
previously, any yawing maneuver will result in an unstable action of the 
booster in the yaw direction. That is, the direction of rotation of the yawing 
moment acting on the booster due to the orbiter exhaust plume impingement 
will be in the same direction as the yawing maneuver. Hence the yawing 
moment due to the orbiter exhaust plume impingement produces an unstable 
rotation about the booster z axis (Fig. 2 -lb). This action could result from 
any one of several causes. These include: (1) one engine out of the (side- 
by-side" engine arrangement; (2) an engine lag in the ignition sequence for 
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the side-by-side mounted engines; (3) angularities and misalignment of the 
wings, horizontal and vertical tail sections mating to the booster fuselage; 
and (4) angularities and misalignment in the engine manufacture and attach- 
ment at the engine gimbal station. 

The second item of consequence is the roll moment acting on the vehi- 
cle during a vehicle rolling or yawing maneuver. During the initial phase of 
the staging sequence a rolling maneuver in the positive roll direction (Fig. 2-lc) 
will produce a moment which tends to increase the roll. As the vertical dis- 
tance as well as the axial distance between the vehicles increases, however, 
the plume external shock will impinge on the v. .ng or horizontal tail surfaces. 
This greatly increases the local loading on the surface which is being rotated 
into the plume. Hence, the increased loading will produce a change in the 
roll moment direction and will tend to then act in a direction opposite the 
vehicle roll direction. It should also be noted that a yawing maneuver will 
produce a roll moment acting on the vehicle. This results since the yawing 
rotation about the vehicle z axis will tend to produce non symmetric or biter 
exhaust plume impingement. The loading on the wings horizontal and vertical 
tails are then unequal so that a roll moment is also produced. 

Vacuum conditions were assumed to exist on all surfaces not subject 
to the orbiter exhaust plume impingement. A preliminary calculation indi- 
cated the aerodynamic force at the separation altitude of 283,000 feet is 
negligible compared to the predicted impingement forces. 

The results of this study were obtained using an equivalent plume which 
has the same mass flow rate as the main stage engines. This analysis included 
the effects of the nozzle lip shock in the impingement analysis. A more de- 
tailed study which includes the effects of the nozzle internal shock wave and 
the p boc k wave resulting from the intersection of the two engine plumes will 
result in a redistribution of the local impingement properties in the vicinity 
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of the shockwave. However, the trends predicted by this study are con- 
sidered to be a reasonable representation of the actual case. Therefore, the 
results of the study during this phase should give realistic data for use in 
preliminary staging studies and for defining potential problems such as those 
discussed above. 

2.2 IMPINGEMENT HEATING ANALYSIS 

Details of the methods used to calculate the plume impingement heating 
rates are presented in Appendix D. Briefly, yawed infinite cylinder theories 
were used to predict the stagnation line heat transfer rates to the fuselage, 
wings and tail surfaces in the continuum regime. Simple heat transfer dis- 
tribution functions which are st ongly dependent on the ratio of the local 
pressure to the stagnation line pressure were used to predict heat transfer 
rates off the stagnation lines. For the range of separation distances studied, 
only the continuum regime was encountered. 

The vertical tail (sketched below) will experience the highest impinge- 
ment heating rates. The point A noted on the sketch is at the leading edge of 
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a chord located 27 inches from the top of the tail. The heating rates to this 
point are presented in Fig. 2-24 as a function of orbiter /booster separation 
distai.ee. Shown on the figure is a line separating the turbulent from the 
laminar heating rates. The heating rates presented in Fig. 2-24 are full 
thrust values. Heating rate values for partial thrust values can be obtained 
from the following equation using the chamber pressure ratio, i.e.. for tur- 
bulent flow. 

^turb " ^full 

thrust 

and for laminar flow, 

^lam °full 

thrust 

Fig. 2-25 shows the variation of the heating rate along the stagnation 
line on the vertical tail for cases*l-14 denoted in Fig. 2- In. For case 1, the 
variation of the heating rate along the chord of the vertical tail containing 
point A (refer to Fig. 2-23, page 2-6) is plotted in Fig. 2-26. It can be seen 
in Fig. 2-26 that the heating rate drops off rapidly away from the stagnation 
line. 


/ Chamber Pressure ) 
V 3000 psi ' 


0.8 


( Chamber Pressure 
' 3000 psi ' 


,0.5 


Shown in Fig. 2-27 are the heating rate contours on the vertical tail for 
cases 1, 4, 10, 12 and 15. The heating rates shown in Fig. 2-27 do not in- 
clude any localized effect due to the nozzle lip shock/fin bow shock inter- 
action phenomena. It is planned to trea these phenomena in more detail in a 
subsequent report. Fig. 2-27e shows the increased heating to the vertical tail 
as the orbiter engine exit plane is moved closer to it. The variation of the 
heating rates along the top of the fuselage ahead of the root of the tail are shown 
in Fig. 2-28. They also do not contain any localized effects due to the nozzle 
lip shock/fuselage bow shock interaction phenomena. The distributions of the heating 


* Case denotes Trajectory Point on Fig 2- In 
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rates away from the stagnation line for a chord o.» the fuselage 5 inches in 
front of the root of the vertical tail are shown in Fig. 2-29 for cases 10 and 
12 . 


In order to obtain the effect of these impingement heating rates on the 
skin temperature o f the vertical tail a transient one -dimensional heat con- 
duction analysis was performed. The separation trajectory used was the one 
labeled "trajectory 1, nominal" in Fig. 2-ln. The thrust -time buildup 
curve used is shown in Fig. 2-30. The points on the vertical tail that were 
analyzed are points A and B noted in Fig. 2-23. At point A, the leading edge, 
a 0.25-inch thick, carbon/carbon material was used. This leading edge 
material extends back to point B (approximately 15% of the chord length) 
where a titanium skin, 0.05-inch thick was used. The heating rates to points 
A and B versus time are plotted in Fig. 2-31 for this trajectory. The surface 
temperature -time histories are shown in Fig. 2-32 for the leading edge (point 
A) and the side of the vertical tail (point B). Point A reached a maximum 
temperature of 1800°F and point B reached c maximum temperature of 480°F, 
both below the structural temperature limits of 2700°F and 600°F, respec- 
tively. An initial temperature of 300°F was assumed in this analysis. 

Since the more vertical separation trajectory maneuver is also of in- 
terest, the separation trajectory labeled "trajectory 2, vertical" in Fig. 2-ln 
was also analyzed. The heating rate contour plots on the vertical tail for 
selected orbiter/booster separation positions along this trajectory are shown 
in Fig. 2-33. 

The variation of the heating rates along the leading edge of the vertical 
tail for selected points along trajectory 2 are shown in Fig. 2-34. The heat- 
ing rate -time histories for points A and B on the vertical tail for trajectory 
2 are shown in Fig. 2-35. The surface temperature -time histories for these 
points are shown in Fig. 2-32. It can be seen in Fig. 2-32 that the leading 
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edge temperature reached the same maximum (1800°F) as in separation 
trajectory 1. The titanium skin {point B) reached a temperature of 440°F 
or 40°F less than trajectory 1. 

For separation trajectory 1 there was very little direct plume impinge- 
ment on the booster fuselage ahead of the root of the vertical tail, therefore 
no surface temperature response was calculated for the booster. Only cases 
10 and 12 had any impingement on the booster and these heating rates were 
rather low (4-6 Btu/ft^ -sec). This, combined with the fact that the heating 
rate pulse would only be for 2 to 3 seconds indicated that the surface tempera- 
ture would remain below the structural limit of 600°F. For separation tra- 
jectory 2 there was no direct plume impingement on the booster ahead of the 
root of the vertical tail. 
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Schematic of Space Shuttle Booster Geometry and Booster 
Representation for Plume Impingement Analysis 
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Note: 1. Drawing not to scale 

2, Force and moment arrows 
indicate positive directions 

3. Booster reference coordinate 
system located at booster nose. 


Fig. 2.1-b - Three-Dimensional Perspective of the Booster Reference 

Coordinate System and the Orbiter Plume Coordinate System 
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Transformation matrix from the booster coordinate system to the 
..roller engine coordinate (order of axis rotations are: v. t, , c?) 
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Fig. 2-lc - Three-DimensionaL Perspective of Rotations About the Booster Axes as 

used in the Booster Pitch, Roil and Yaw Plume Impingement Calculations 
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Fig. 2 -lg — Three -Dimensional Perspective of th>~ Booster Right Wing 







2-18 


LOCKHEED • HUNTSVILLE RESEARCH & ENGINEERING CENTER 

w 

n 



-lh — Three-Dimensional Perspective of the Booste r Left Wi 




i 

I 

* 


4 

2-19 I 

LOCKHEED • HUNTSVILLE RESEARCH A ENGINEERING CENTER 


Fig. 2 -li — Three-Dimensional Perspective of the Booster Left 
Horizontal Tail 
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Fig. 2 - 1 j ~ Three-Dimensional Perspective of the Booster 
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Fig. 2-lk — Three -Diminsional Perspective of the Booster Vertical Tail 
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Fig. 2.1-n - Graphic Representation of Data Points Used to Calculate the Booster 

Plume Impingement Data with Separation Distance Between the Orbiter 
and Booster Vehicles 
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Fig. 2-4 — Variation of the Booster Moment due to the Orbitcr Exhaust 

Plume Impingemen* with Orbiter Separation from the Booster 
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Fig. 2-7 — Variation of the Booster N-. lent, due to the 

Orbiter Exhaust Plume Impingement with Orbiter 
Separation Distance from the Booster 
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Fig. 2-8 Variation ox the Booster Plume Impingement Force with 
Orbiter Separation Distance from the Booster 
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Fig. 2-9 — Variation of the Booster Plume Impingement Force with 
Orbiter Separation Distance from the Booster 
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Fig. 2-11 — Variation of the Booster Plume Impingement Force with 
Orbiter Separation Distance from the Booster 
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Fig. 2-12 — Variation of the Booster Plume Impingement Force with 
Separation Distance from the Booster 















LMSC/HREC D162657 





Fig. 2-13 — Variation in the Booster Force due to Orbitei 

Exhaust Plume Impingement with Orbiter Separation 
Distance from the Booster 
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Fig. 2-15 — Variat'on of the Booster Moment due to Plume Impingement 
with Separation Distance of the Orbiter from the Booster 
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Variation of Booster Roll Moment Due to Orbit Exhaiim Plume 
Impingement w’th Orbiter Separation Distance from tin Booster 
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-24 - Variation of the Impingement Heating Rate to the Booster Vertical 
Tail at Point A as a Function of Orbiter Engine Exit > lane 

Booster Separation Distance 









Fig. Z-Z 5a - Variation of the Heating Rate Along the Stagnation Line on the Rooster Vertical Tail 












Fig. 2-26 — Variation of Heating Rate Along Chord on Vertical Tail for Case l(Fig. 2-ln) 
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Fig. 2-27d — Heating Rate Contours on the Vertical Tail for Case 12 
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Fig. 2-29 — Variation of Fuselage Heating Rates off the Stagnation Line at a Point 
Immediately in Front of the Root of the Vertical Tail. 
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Fig. 2-33d — Heating Rate Contour Plats for the Vertical Tail for Case 19 (Fig. 2- In) 
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Fig. 2-33e — Heating Rate Contour Plats for the Vertical Tail for Case 20 (Fig. 2- In) 
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Fig. 2-34 — Variation of the Stagnation Line Heating Rate on 
the Vertical Tail foi Cases 16. 17, 18 and 21 
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Fig. 2-35 — Heating Rat^i, Time History for the Vertical Tail 
for Separation Trajectory No. 2 
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Section 3 
CONCLUSIONS 

In this study, the plume impingement heating rates, surface tempera- 
tures, forces and moments incurred by the MSC booster during the staging 
maneuvers were calculated. Although this study was preliminary, the data 
presented in this report can be used to determine a staging separation trajec- 
tory which will result in acceptable impingement surface temperatures, 
forces and moments. The data include the effects of roll, pitch and yaw 
maneuvers for the booster and are presented in graphical form. 

This study is considered preliminary because: (1) the basic exhaust 
plume used to simulate the multiple orbiter engine exhaust plume was preli- 
minary in nature and reflects the simplifications and assumptions used to 
generate it, and (2) the localized effects of plume shock/bow shock inter- 
action phenomena on heating rates were net trea> ^d. 

Due to the multitude of plume shock/bow shock and fuselage bow shock/ 
vertical tail bow shock interaction situations, it is advi able to first select a 
staging separation trajectory neglecting these interactions, and then look at 
these interactions in more detail to determine the optimum staging separation 
trajectory. Each interaction is a unique situation and must be analyzed 
separately. 

In order to determine if a staging separation trajectory resulted in 
"acceptable" heating rates, it was necessary to perform a transient heat con- 
duction ana’vsis to determine if the booster skin temperature exceeded the 
maximum structural temperature limits. Using a simplified thrust buildup 
curve, the heating rate-time histories to selected points on the MSC vertical 
tail were used in a simple one -dimensional heat conduction program to cal- 
culate the surface temperature -time histories. These surface temperature - 
time histories were calculated for two separation trajectories. The maximum 
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calculated surface temperature at the leading edge of the vertical tail was 
1800°F for both of the separation trajectories analyzed. The titanium skin 
on the sides of the vertical tail reached 480°F and 440°F for the two separa- 
tion trajectories analyzed. The temperatures are beluw the maximum 
structural temperatures of 2700°F for the NASA-MSC carbon/carbon leading 
edge material and 600°F for the tita lam skin on the side of the vertical tail. 

The predicted forces and moments resulting from the orbiter main 
stage engine exhaust plume impingement on the booster give realistic data 
for use in preliminary studies concerning the staging sequence of the space 
shuttle vehicle. These data, though of a preliminary nature, are considered 
to predict closely the force and moment trends to be experienced by the 
booster due to orbiter exhaust plume impingerr. nt. A detailed analysis which 
includes the effects of the orbiter nozzle internal shock wave and the shock 
waves produced by the two engine plume interaction region will tend to re- 
distribute the impingement flow properties locally. However, the integrated 
force values should remain substantially unchanged. 

The present study disclosed several areas which could pose potential 
problems as a result of the orbiter exhaust engine plume impinging on the 
booster. These are: 

1. A positive pitch rotation (rotation about the booster y oxis; Fig. 
2-lc) will tend to produce a booster pitching moment which acts 

in the same direction as the booster rotation if the plume impinge- 
ment center of pressure is aft of the booster c. g. This could be 
somewhat critical during the initial phases of the staging sequence 
since in this case the impingement will always be aft of the booster 
c. g. 

2. A yawing rotation of the booster about the Z axis (Fig. 3-lc) will 
produce an unstable effect in that the moment which will act on the 
booster in the yaw direction will augment the booster yaw motion. 


3-2 


LOCKHEED ■ HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC/HREC D162t>57 


Hence the booster would continue to rotate since the large vertical 
tail would become more fully exposed to the exhaust plume. A yaw- 
ing rotation will also produce a roll moment acting on the vehicle 
which will be acting in the direction of vehicle roll. Hence, if the 
moments resulting from plume impingement were to become sub- 
stantia), the vehicle would tend to tumble. 

3. Tne deflection of the wing horizontal tail surface as a result of the 
plume impingement load dist; button would result in larger loadings 
and hence increase tne magnitude of the potential problem discussed 
above. 

Areas which should be investigated *n addition to the present analysis are 

(1) the engire-out case for side-by-side mounted engines (since this would pro- 
duce moments which wil 1 tend to rotate the vehicle about the roll and yaw a>es) 

(2) effects of misalignment of the wings and tail surfaces with respect to the 
booster fuselc >; and (3) misalignment of the orbiter main stage engine with 
respect to the gimbal plane. 
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Appendix A 

Representative plumes for the single and multiple nozzle Space Shuttle 
orbiter engines were generated for the analysis of Phase I of this study. 

The plume analysis considered the Rocketdyne version for the orbiter main 
stage engine (Table A.l). 

A. 1 THERMOCHEMISTRY ANALYSIS 

The proposed Space Shuttle main propulsion engines employ the concept 
of staged combustion which is reflected in the gas thermochemical data 
(Table A.2). Thermochemical properties of the preburner combustion pro- 
cess were generated using the NASA-Lewis Chemical Equilibrium Composi- 

4 = 

tion (CEC) Computer Program. The propellant used in the preborner 
analysis was: as the oxidizer and as the fuel with a mixture ratio of 

6. An estimate of the change in total enthalphy of the preburner combustion 
gases across tl hydrogen and oxygen turbines was obtained from turbine 
power and efficiency data. Thermochemical properties resulting from the 
main burner combustion process were then calculated using the preburner 
combustion products as the propellant and including the adjusted total enthalpy 
in the analysis. The data from this analysis was used in all the subsequent 
flowfield analyses and are presented in Table A.2. 

A.2 PLUME FLOWFIELD ANALYSIS 

Two plume flow fields for the orbiter engines were generated for use in 
the preliminary analysis. These were: (l)the engine operating at 283.000 feet 


4 < 

McBride, Bonnie, and S. Gordon, "Preliminary Description of CEC, A Com- 
puter Program for Calculation of Chemical Equilibrium Compositions with 
Applications," NASA-Lewis Research Center, Cleveland, Ohio, May 1961. 
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with the vehicle Mach number effect on the plume shape considered, and (2) 
the engine exhausting into quiescent conditions at 283,000 feet. The latter 
plume (which is considerably larger in size than the former) will be used to 
establish an upper bound on the booster impingement data and will be included 
in the final report of the subject contract. 

All plume calculations were performed assuming the engine to be operating 
at full t* rust (engine chamber pressure = 3000 psia). The coordinates r,x 
were made non-dimensional with respect to the engine exit radius. Thus, the 
single and equivalent multiple nozzle flow fields are obtained by multiplying 
the field coordinates by the appropriate engine exit radius. 

The plumes were generated using two method- of- characteristics com- 

puter programs. The calculations were initiated at the nozzle exit plume 

with a start line based on a one -dimensional nozzle solution. Real gas effects 

were included in the nozzle solution with the flow angle distributed across the 

exit plane. Input data included: (1) nozzle area ratio of 200/1; (2) lit ingle = 

6.6 degrees and (3) the nozzle external flow condition M = 12.08 at ar altitude 
of 283,000 feet. 

Assumptions used in the calculations are: (1) the flow is inviscid and 
axially symmetric; (2) the flow is rotational; (3) the oxidizer-to-fuel ratio 
(O/F) is constant; and (4) the flow is non-continuum for Knudsen numbers 
greater than "1 /collision number." A lip shockwave (See Figs. A- 1 through 
A- 12) was considered in the analysis. This shockwave is the result 
of the nozzle flow being underexpanded at the exit plane. The shockwave 
follows the plume boundary over the region of interest. 


$ 

Prozan, R. J., "Development of a Method-of-Characteristics Solution for 
the Supersonic Flow of an Ideal, Frozen or Equilibrium Reacting Gas/Mixture, 
LMSC/HREC D1 62220-HI, May 1970. 

Sniit.*, S. D., and A. W. Ratliff, "User's Manual — Variable o/F Method-of- 
Characteristics Program for Nozzles and Plume Flow Analysis," LMSC/ 
HREC D1 62220-1, May 1970. 
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At the nozzle exit plane, the primary constituents of the flow are H^O 
(gas) and H^tgas). The mixture percent composition is 75.6% H^O and 
24.4% H-, (global molecular weight of 14.112 lb m /lb m -mole). Based on a con- 
sideration of the primary constituents, H^O and H-,. The flow approaches a 
chemically frozen condition at a pressure of 0.4 atmospheres which corre- 
sponds to a Mach number of 4.2. For the engine used in the current study 
this condition would occur inside the nozzle. After the flow approaches the 
chemically frozen condition, the isentropic exponent, gamma, then becomes 
a function only of temperature. Because of the lack of thermochemical data 
at low temperatures, those regions of the flow field having Mach numbers 
greater than 3.392 were treated as an ideal gas at a molecular weight of 
14.112 lb m /lb m *mole and an isentropic exponent of 1.3347. 

Non-gaseous forms of the mixture constituents were not considered in 
the thermochemical calculations. The ability of the NASA-Lewis (CEC) 
Thermochemical program to handle liquids and solids is restricted by the 
lack of thermochemical data for non-gaseous species. Because condensed 
species might occur in the plume, vapor pressures for the composition were 
examined to determine where incipient condensation might exist. Conditions 
are not favorable for the condensation of I^O or 1^ over the range of interest 
for the current problem. This is due to the high combustion pressure which 
results in a high static pressure in the plume which delays the on-set of 
condensation. 

The flow fields for the respective plumes are described by the streamline 
distribution shown in Figs. A-l and A-7. The mass flow associated with each 
streamline is indicated with the mass flow given as a percentage of the total 
mass flow in the plume. The plume flow fields are mapped by contour plots 
of Mach number, local static pressure, temperature, density, and pitot total 
pressure as shown on Figs. A-2 through A-6 and A-8 through 12, respectively. 

Knudsen number calculations were made at various points in the plume 
flow field co determine if the flow were in the free molecular flow regime. 
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Temperature gradients indicate that the flow is well within the continuum 
flow regime. 

To consider the effects of the engine ignition transients the local flow 
properties should be adjusted by the ratio of the chamber pressure at time, 
t, in the ignition sequence to the chamber pressure of 3000 psia at full thrust 
conditions. 
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Table A-l 

OPERATING CHARACTERISTICS FOR THE ORBITEL 
MAIN STAGE ENGINE 


Engine Manufacturer: Rocketdyne 


Preburner 


Oxidizer: 0 7 
Fuel: H 2 


Mainburner 


Oxidizer: O^ 
Fuel: H 2 , H-,0 


Engine Parameters 

A/ A* = 200/1 
O/F =6.0 

Effective equivalence ratio = 1.3227 
(main burner) 

P^ = 3000 psia 

0 = 6.6 degrees 

exit 
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Table A-2 

THERMOCHEMICAL DATA USED IN THE FLOW FIELD ANALYSIS 
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Fig. A-l - Distribution of Mass Flow Rate in the Orbiter Main Stage Engine 
Plume with Hypersonic Plume External Flow Conditions 
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Hypersonic Plume External Flow Conditions 
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Fig. A-6 - Pitot Total Pressure Contour Map of the Orbi+er Main Stage Engine 
Exhaust Plume with Hypersonic Plume External Flo*v Conditions 















BStE 

LLuV 


o a 

So «•• 

;?■ 3 - 
• 

v 'O • m 

* § I 

* -a | • S o 

S-2 S ; . 
i • a o.S E 

• *• 3 * 3 

; gi • " * 
i. u E« 3 " 

: 55 8 5 1 

i HOZ« N 


; 



1 






V 






J 

% 





1 

A] 






’iu 






1 

■ 




- 


\ 

'[ 


























































■ J 






















i 

■ 

- 


— 























X-U-Ut 



Fig. A -8 - Mach Number Contour Plot for the Orbiter High Altitude Plume 
Fxhausting into Quiescent Conditions 
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Fig. A -9 - Static Pressure Contour Map for a High-Altitude Orbiter Plume 
Exhausting into Quiescent Conditions 




Fig. A- 10 - Static Temperature Contour Map for the Orbiter High Altitude 
Plume Exhausting into Quiescent Conditions 
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Symbols 


c 

dL 

f 

i, j, k 
n 


/S /\ A* 

nx, ny, nz 
r 

s 


S 

t 

r 


x, y, z 
x' , y' , z' 


NOMENCLATURE 


chord , in. 

differential length vector 
airfoil cross section function 
unit vectors subshape system 

unit normal 

defined in text 
root chord, in. 
span , in. 
surface function 

taper ratio (wing tip chord/root chord) 

coordinates of point in local subshape system in. 
local leading edge coordinate system (normalized) 


Greek 

a 

3 

0 


angle between normal to surface and y coordinate 

defined in text 

sweep of y' = constant line 

leading edge sweep angle 


i 
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Appendix B 


INTRODUCTION 

A comprehensive plume impingement program has been developed 
under several support contracts. The impingement analyses conducted to 
date were mainly concerned with orbiter vehicles and the library of coded 
subshapes (conic sections, rectangular, triangular and circular flat plates) 
were adequate to describe the vehicle's geometry. Under Contract NAS9- 
11156, however, it was necessary to describe airfoil subshapes in order to 
accurately simulate the Space Shuttle vehicle geometry. 

This Appendix describes the theoretical approach to this airfoil descrip- 
tion. The airfoil considered for this study was NACA 0012-64 which is sym- 
metrical. The subshape calculation has been coded, however, to permit the 
use of any airfoil as long as there is no spanwise variation in the normalized 
contour. 

DISCUSSION 

Consider the airfoil shown in Fig. B-l. Let x' = spanwise location 
normalized by the span; y' = chordwise distance aft of leading edge norm- * 

alized by the chord; z' = airfoil thickness normalized by the chord. The 
airfoil cross section equation is 

S (x 1 , y' , z') = z' - f(y') = 0 


* 

Wojciechowski, C. J. and A. W. Ratliff, "Study of High Altitude Plume 
Impingement," LMSC/HREC D149118, August 1969. 
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from which i'« follows that 


VS = - 


9f_ 

9y' 


j + k 


It can be seen however that 


-Q-p = tan (a - y) = - cota . 


Consider the change of variables 


X = s x' ; y = X tanQ + r 1 - (1 - t ) 
s r 


1 - y'; z=r l-(l-t r )i z' 


where 


c = r 


[i-u-t r )f] 


Now 



= tan 0 = tan 0 
y' = const s 


j: 

s 


( i - V y 


and 


dz 

dx 


y* = const 


- f (1 -t r )z' = tan 3 . 


From the chain rule 



as 8 H 
^ si- 


where 


Bx 1 _ _ 1 _ 9y< tanQ . 9z' tanB 

9 x ~ s ’ 9x c ’ 9x c 
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9 x 1 
9y 



_1_ 8 z 

c : 8 


9x' 

9z 


0; 


I? •« 


8 z 1 

9 z 


The partial derivatives of S with respect to are therefore; 

= _ cota tan0 . tang 

9x 

8_S 
9y 

8S 
8 z 

The normal to the surface is 

— ^ — - 
n=n i + n j + n k 
x y J z 

2 ~ 2 ~ 2 
w n + n + n 
V x y z 

where 

n = - (cosa tan© + sina tan p) 

n = cosa 

y 

n = sina • 
z 


c c 

cotg 

c 

J_ 

- • 

c 
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1 

Elemental Area 


1 

Taking the differential of the transformation yields 


» 

■ 

dx = s dx' ; dy = s(tan0+ t?: " 1 dx 1 + c dy' ; dz = s tan|3 dx' + cdz' . 



So that a differential length vector may be written as 



dL = sdx* i ++ |s tan0 dx* + c dy'j j + (s tan 3 dx' + cdz') k . 



If a step is taken in dx 1 (dy 1 = 0 ) 



dLj = s £i + tan0 j + tan(3 k J dx' + cdz' k . 


1 

If a step is taken in dy' (dx* = 0 ) 


1 

dL <2 = c iy' j + cdz 1 k . 


■ 

1 

To ensure that these differential steps are taken along the surface - .t is 
necessary that 



*i • dLj = n • dL^ = 0 . 


1 

• 

But 


\ n 

n *dL, = dz 1 = 0 


1 

so that 


1 

dLj = s £i + tan0 j + tanp k J dx' . 


1 

\ 
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It was found that the calculation of dz 1 from n 'dL^ suffered from 
numerical inaccuracies. The quantity dz* is therefore found directly from 
the function and 


d l 2 = c 0 + 371- k ) d y’ • 


The differential area is 


whence 


dA = dL 1 x dL 2 


s c 


4^dy' = [i + tantj + tanp kj x 


T + — k 
dy* 


dA 


s c dx 1 dy 


• = [ tan 6 Ip * tan p] i ” fp J + k 


and finally 


dA = 


tan0 If! ‘ tan p| 2 + 


dz 1 

dy’ 


2 1 1/2 

+ 1 sc dx* dv* 


/ 
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Appendix C 

The staging sequence of the Space Shuttle vehicle will result in those 
exhaust plumes of the orbiter main stage engines impinging on the booster 
surfaces. Calculations were made to determine the forces and the resulting 
moments imposed on the booster due to the orbiter exhaust plume environment. 

The impingement analysis mapped the force and moment distribution on 
an area to sufficiently analyze the expected separation trajectories. The data 
include the effect of roll, pitch and yaw maneuvers for the booster and is pre- 
sented in graphical form. 

C.l PROBLEM DESCRIPTION 

The geometrical orientation for the Shuttle arrangement is shown in 
Fig. C.l. A reference coordinate system is located at the nose of the booster 
on the vehicle centerline. All forces and moments presented are referenced 
to this point. A simulated multiple nozzle was used to represent the orbiter 
main stage engines. The simulated nozzle has the same mass flow rate as 
the orbiter engine cluster and approximately the same plume shape. 

In the plume impingement analysis, the various components of the booster 
were represented by four (4) basic subshapes (Fig. C.2). These were: (1) a 
cylinder for the upper portion of the fuselage; (2) flat plates for the lower sides 
of the fuselage; (3) a conical nose and airfoil shape NACA 0012-64 for the wing, 
vertical and horizontal tails. 
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C.2 DEFINITION OF THE PLUME IMPINGEMENT CHARACTERISTICS 

The impingement characteristics were calculated by Lockheed's Plume 
Impingement Computer Program*. Basically, the procedure for obtaining the 
impingement characteristics is as follows. The vehicle composite shape and 
the particular motor are located with respect to the reference coordinate 
system. The composite shape is then represented by a series of geometrically 
simple subshapes (each with its own local coordinate system). These subshapes 
are further divided into smaller incremental areas. Local flow condi- 
tions on a given body are obtained by locating the centroid of each elemental 
area in the free plume** from a search of the plume flow field which is stored 
on magnetic tape. The flow properties thus obtained are assumed to act uni- 
formly over the elemental area. 

The body local flow regimes were determined by the value of the local 
Knudsen number. In this analysis, continuum flow was assumed to exist for 
Knudsen numbers less than 0.01; transitional flow exists for Knudsen numbers 
between 0,01 and 10.0; and Knudsen numbers greater than 10.0 define the 
"free molecular flow" regime. The characteristics lengths used in the 
Knudsen number calculations were a representative dimension (element area 
dimension) of the body in question so that the elemental area local flow regime 
is determined by comparing the mean free molecular path of the local plume 
to the body dimension. The analysis assumed the surfaces to be of an infini- 
tesimally small thickness. 


Wojciechowski, C. J., and A. W. Ratliff, "Study of High- Altitude Impinge- 
ment," LMSC/HREC D149118, August 1969. 

Prozan, R. J., "Development of a Method-of-Characteristics Solution for 
the Supersonic Flow of an Ideal, Frozen or Equilibrium Reacting Gas 
Mixture," LMSC/HREC D162Z20-UI, May 1970. 

Smith, S. D., and A. W. Ratliff, "User's Manual- Variable O/f Method-of- 
Characteristics Program for Nozzles and Plume Flow Analysis," LMSC/ 
HREC D162220-I, May 1970. 
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In the numerical integration scheme, the cylindrical portion of the fuse- 
lage and the conical nose section were divided in rings of 50.0 inches in incre- 
ment with each ring further subdivided into 10.0-degree radial incremt nts. 

All wing sections were divided into approximately 40.0 inch increments in the 
span direction. As discussed in Appendix B, a cosine distribution function is 
used to distribute the data points in the chord direction. This scheme biases 
the point distribution toward the wing -leading edge where largest grad: nts in 
surface curvature occur. Twenty (20) increments in the chord di: ac'-i n were 
used on the wing and horizontal tail surfaces while thirty (30) inci era its were 
used on the vertical tail surfaces where the highest values of hef ting rates and 
impingement force distribution occur. Uniform flow was assumed to exist over 
each elemental area of the respective surfaces with vacuum conditions considered 
to exist on the surfaces not exposed to the engine exhaust plume. Cha acteristics 
dimensions used in the Knudsen number calculations to determine tht urface 
local flow regime were the element area dimensions and the diameter ■■■ and 

lengths of the conical shaped bodies. Also, a representative plume molecular 

-6 

diameter of .1 x 10 in. was used in the Knudsen number calculations. It 
should be noted that the force integration results reflect the assumptions used 
to generate the engine exhaust plume. * 

C.3 PLUME IMPINGEMENT FLOW REGIMES 

In the continuum regime, the local impact pressures were cal :ulated 
from a modified Newtonian impact theory as given by 

P. =P (1 + £ KM 2 sin 2 ct) (C.l) 

imp oo ' 2 n 

with K given by 

k = c p (.8i4+-M|-\ (c.r- 

s' a / 

*Prozan, and Smith, LMSC/HREC D162220-I and LMSC/HREC D1622<i0 -.III, 

May 1970. 

Buck and McLaughlin, ASP-TDR-63-552, August 1963 (Confidential). 
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Fundamental Newtonian theory assumes that the force felt by the body 
is due to the time rate of change of the normal momentum and is independent 
of the angle of incidence of the body surface so that there is no transfer of 

the tangential momentum. The parameter K is included to account for the 

$ 

higher pressures noted experimentally for small -incidence angles on flat 
plates. For this analysis, K ranged from 1.8 to 6.0 and was set to 6.0 for 
local incidence angles less than 3 deg. The forces and moments were then 
obtained by numerically integrating the local impact pressure distribution 
over the body. 


The free molecular flow regime is characterized by the fact that the 

local properties on the body depend only upon collisions of the gas molecules 

with the body surface. A Lockheed/Huntsville- developed impingement 
** 

program assumes environmental equilibrium so that a Maxwellian- type 
flow exists and the local flow properties are calculated from kinetic theory 
and Maxwell's law of distribution of molecular velocities. The local force 
on the elemental area is obtained by computing the forces produced by the 
incident and reflected molecules and adding them. The force is given by 


dF 

dA 2(3 



(ke + ly + tq) 


+ y (1 + erf (yS) + j- 





y/T( 1 + erf (yS) + Se" 7 



2 ,. 2 "| 

V£“{1 + erf(yS) + e~ 7 H 


(C. 


Buck, M. L. and E. J. McLaughlin, "A Technique for Predicting Pressure 
Distributions and Aerodynamic Force Coefficients for Hypersonic Wing 
Reentry Vehicles (U), n ASP-TDR-63- 552, Wright-Patterson AFB, Ohio, 
August 1963 (Confidential). 

$ $ 

Wojciechowski, C. J., and A. W. Ratliff, "A Study of High Altitude Impinge- 
ment," LMSC/HREC D149118, August 1969. 


C-4 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 


LMSC/HREC D1 62657 


The total force in the free molecular regime is obtained by numerically 
integrating Eq. (C.3) over the body surface. A detailed development of Eq. 
(C.3) is given in a Lockheed/ Sunny vale document? 

In the transition flow regime the local force is calculated from the 
following empirical expression from a Lockheed/Huntsville document.* 

P tran = P co.t + < P f.„,. ' P co„,> 8i ” 2 [' l 1 / 3 + lo 8l0 K „] • < c ' 4 > 

Equation (C.4) is a function of Mach number and Knudsen number and yields 
a smooth transition from the continuum to the free molecular flow regime. 

C.4 INTEGRATED PRESSURE FORCES AND RESULTING IMPINGEMENT 
MOMENTS 

The integrated forces for the booster are obtained by numerically inte- 
grating impact pressure distributions over the individual booster subshapes. 
The force on an elemental area of a subshape is given by 


F f= P imp dAl l* < C ‘ 5 > 

The resultant force acting the booster due to orbiter exhaust plume impinge- 
ment is then obtained by summing the forces over each subshape. The 


Sentman, Lee H., "Free Molecular Flow Theory and Its Application to the 
Determination of Aerodynamic Forces," LMSC-448514, Lockheed Missiles 
& Space Company, Sunnyvale, Calif., October 1961. 

Ratliff, A. W., et al., "Analysis of Heating Rates and Forces on Bodies 
Subject to Rocket Exhaust Plume Impingement," LMSC/HREC A791230, 
March 1968. 
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resultant force acting on the booster and referenced to the booster nose is 
then expressed as 




(C.6) 


Since the flow distribution is assumed to be uniform over each elemental 
area of a given subshape, the resultant force for the elemental area chen acts 
at the centroid of the element. The moment acting on the element area is then 
given by the expression 


M i = r f xF i , (C.7) 

where the vector is the distance vector from the origin of the subshape 
coordinate system to the centroid of the element in question. The resultant 
moment acting on the booster is given by a summation of the moment distri- 
bution over the booster. The mathematical expression becomes 




(C.8) 


The impingement calculations assume the flow field to be inviscid. There- 
fore any boundary layer effects which produce plum • impingement outside the 
calculated plume boundary are disregarded as is any vise -us (skin) forces 
between the plume and the respective subshapes. 
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Plum's 

Boundary 


.Shock Wave 
Or biter Engine 



Drawing not to scale 

Force and moment arrows 
indicate positive directions 

Booster reference coordinate 
system located at booster nose. 


Fig. C-l " Three-Dimensional Perspective of the Booster Reference 

Coordinate System and the Orbiter Plume Coordinate System 
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Symbols 

°P 

g 

H 

h 

k 

M 

P 

Pr 

» 

q 

R 

T 

t 

u 

V 

v 

X 

Greek 

0 

y 

6 


NOMENCLATURE 

specific heat (Btu/lbm-°R) 
gravitational constant (ft- lbm/lbf-sec 2 ) 
total enthalpy (Btu/lbm) 
static enthalpy (Btu/lbm) 
thermal conductivity (Btu/Ft- sec~°R) 

Mach number 
pressure (lbf/Ft 2 ) 

Prantl number 

convective heating rate (Btu/Ft^"sec) 
radius of cylinder (Ft) 
temperature (°R) 
time (sec) 

normal component of freestream velocity (ft/ sec) 
velocity (ft /sec) 

tangential component of freestream velocity (ft/ sec) 
length (ft) 

angle of attack - radians 

ratio of specific heats 

flow deflection angle - radians 
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P density (slugs/Ft^) 

H viscosity (slugs/Ft- sec) 

Subpcripts 

e boundary layer edge values 

o frustum stagnation conditions 

r recovery value 

s stagnation line conditions 

w wall conditions 

00 free stream condition ahead o* shock 

1 property upstream of shock 

2 property downstream of shock 
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Appendix D 

This appendix discusses the heac transfer theories used in this study 
and contained in the plume impingement program (Ref. D. 1). For the range 
of separation distances considered in this study only continuum regime heat 
transfer theories were used. The boundary layer convective heat transfer 
theories are discussed in terms of tv o regions; i.e., the stagnation regions, 
and the non-st' -ion regions. 

In the coi -n regime, the MSC booster will experience turbulent 

as well as laminar boundary layer convective heating. For close orbiter/ 
booster separation distances, the local freestream Reynolds numbers foi the 
impinging flow are high and therefore turbulent boundary layer heat transfer 
theory was used. For large orbiter/booster separation distances, the local 
freestream Reynolds numbers of the impinging flow are low, and therefore, 
laminar boundary layer theory was used. The criteria used to establish the 
state of the boundary layer, i.e., either laminar or turbulent, are also pre- 
sented in this Appendix. 

The following sketch denotes the orbiter engine impingement flow field 
on the booster. Noted on the sketch are the flow stagnation lines on the lead- 
ing edges of the wing and tail surfa^** and top of the fuselage. Along these 
lines, stagnation line heating rate methods were used. For heating rate 
distributions off the stagnation lines only the component of the flow in the 
crossflow direction was used. This resulted in simple ratios of iocal-to- stag- 
nation line heating rates which were strongly dependent on the ratio of local- 
to-stagnation line pressures. Noted in the sketch are several possible shock/ 
shock interaction areas caused by the fuselage shock envelope intersecting 
the leading edg’ shock envelopes. 
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Fuselage Impingement 
Shock/Leading Edge Shock 



The effects of these shock interactions were treated ir. this prelimi’ ary 
study in the context ox promoting early boundary layer transition. The in- 
creased heat transfer in these areas due to vortex impingement will be stu. ed 
in a subsequent report. 

Preceding the discussion on the boundary layer convective heat transfer 
theories will be a section on the methods used to obtain the local impingement 
flowfield properties which are used as the inviscid conditions at the edg<=- of 
the boundary layer i*. the heat transfer equations. 

Following the discussion of the boixndary layer transition criteria will 
be a section discussing the methods used to obtain the structural thermal 
response to the impingement heating rates. 
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D.l LOCAL FLOWFIELD CALCULATION METHOD 

This section presents the method used to calculate the local flowfield 
properties along the stagnation lines on the booster fuselage, wing and tail 
surfaces. Inviscid analysis methods were used and no boundary layer inter- 
actions were treated. 

The impingement pressure on the body was calculated using Modified 
Newtonian theory as discussed in Appendix C. With the pressure along the 
stagnation line calculated in this fashion, the remaining local flow properties 
at the boundary layer edge are computed as the real gas conditions behind 
an oblique shock wave which is parallel to the local body slope. 

Basically the solution involves the solution of the continuity, momentum, 
and energy equations plus an equation of state. The following sketch shows 
the flowfield geometry. 



V l =V 2 =V w (D.l) 
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The pertinent equations are 







Conservation 
of Momentum (D.2) 





1 - 



Conservation 
of Ene rgy 


(D.3) 


p l U 1 = P 2 u 2 


Continuity (D.4) 


tan (8-6) ^2 

tan 8 U j 


Oblique 

Shockwave 

Geometry (D.5) 


whe re M , = hi sin 0 . 

1 oo 

An iteration procedure is used to obtain the solution to the preceding 
set of equations. After the shock solution is obtained, the flow is isentrop- 
ically decelerated such that equals zero with a corresponding increase 


D.2 BOUNDARY LAYER HEAT TRANSFER THEORY 


In this section, the boundary layer heat transfer theories are discussed 
in terms of two regions; i.e., the stagnation regions, and the non -stagnation 
regions. Both turbulent and laminar boundary layer theories are discussed. 
The heating rates along the stagnation lines are calculated using yawed in- 
finite cylinder theory. 

D.2.1 Stagnation Regions 

Turbulent Boundary Layers. The heat transfer rates along the stagna- 
tion lines for turbulent boundary layers were calculated using a modified 
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version of thr' turbulent swept cylinder theory of Beckwith and Gallagher 
(Ref. D. 2); i.o.. 


0.0288 

Pr 

w 



(p n ) 0 8 (H - h ) 
w w r w 


(D.6) 


Modifications to this theory include the introduction of the enthalpy 
driving potential to make it more applicable to boundary layers in plumes. 

The complete derivation of Eq. (D.6) is given in (Ref. D.2). Briefly, the 
integral energy and momentum equations for the turbulent boundary layer on 
a yawed infinite cylinder are solved using simplifying assumptions for skin 
friction, velocity and thermal profiles, secondary flow and Reynolds analogy. 
Eq. (D.6) is based on the assumption of a Blasius flat plate skin -friction law 
for the flow along the leading edge and zero secondary flow which is valid 
only at the stagnation line. A comparison of the results of this modified 
theory with test data is given in Ref. D.3. A typical figure from this work 
is shown in Fig. D. 1. 


The transport properties used in the heating rate equation are calcula- 
ted for the individual gas species using the Lennard -Jones potential model. 
The mixing laws for obtaining the viscosity and conductivity of the gas mix- 
ture are the same as those described in Refs. D.5 and D.6. The velocity 
gradient can be estimated from Newtonian theory or from the calculated 
pressure gradients in the stagnation region (Ref. D.4). 

If the crossflow velocity component is supersonic, then the velocity 
gradient is computed by combining the edge of the boundary layer momentum 
equation with the Newtonian pressure distribution for a cylinder, i.e., 



(D.7) 
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If tne crossflow velocity component is subsonic, then the velocity 
gradient is computed by the following equ \tion. 


dV 

dR 


V sin 0 
ao 

2 R 


5.1 cos sin 0) - 1.6 . 
oo 


(D. 8) 


Laminar Boundary Layer The heat .rLnsfer rates along the stagnation lines 
for laminar boundary layers were calculated using the empirical correla- 
tions of Ref. D.7 . The resulting heat transfer correlation obtained was 


4 


0.565 

pT^ 

w 



(D.^) 


This equation agrees closely with the theoretically derived heat transfer 
equation of Beckwith (Ref. D.26). 

D.2.2 Non -Stagnation Regions 

Equations (D.6) and (D.9) are valid only along the flow stagnation lines. 
To obtain heat transfer distributions off the stagnation lines use is made of 
empirical correlations. 


Turbulent Flow. The heating rates off the stagnation lines for turbulent flow 

were calculated using a modified form of the heat transfer ratio derived in 
Ref. D. 2, i. e. , 

~l 0.375 


■ft) 


2ZzL 
2 Y 


&) 

21 iL 1 

y Ts 

' To 


1 - Zl 

To J 


(D. 10) 


Examination of turbulent heating rate distribution data for wed cylinders 
revealed that equation D. 10 correlated the data better than me original form 
of the equation presented in Ref. D. 2. 
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Laminar Flow. The heating rate at a point on the surface for an isothermal 
wall with local similarity assumed can be expressed in terms of the heat 
transfer rate at the stagnation line by, (Refs. D.8 and D.9). 



(D- ID 


The above equation consists of the product of inviscid and viscous quanti- 
ties, The inviscid quantities require a knowledge of the local surface pressure 
distribution as well as the thermodynamics of the inviscid flow. The viscous 
quantities, i.e., g' (0) and g' Q (0), require the appropriate solution to the 
boundary-layer equations which are implicitly related to the inviscid flow 
through the local surface pressure and pressure gradient. The preceding 
statements infer a strong dependence of heating rate ratio Eq. (D ll) on the 
pressure ratio. With this thought in mind, the authors of Ref. D.8 plotted heat- 
ing rate ratios q/q , versus pressure ratio P/P . The resulting plots repre- 
sented an equation of the form 



The data indicate that on hemispherical forebodies the exponent m in 
Eq. (D. 12) assumes values between 0.8 and 0.9 while over the afterbody, m 
assumes values between 0.9 and 1.0. These results demonstrated the strong 
dependence of heating rate ratio on pressure ratio and are useful for obtain- 
ing reasonable estimates of the local heating rates once the stagnation heat- 
ing rates are known. Examination of heating rate ratios for yawed cylinders 
revealed the same dependence of heating rate ratio on pressure ratio. The 
value of m in Eq. (D. 12) which correlated the heating rate ratio best was a 
value of 0.85. Fig. D.2 shows two such correlations over cylinders at angles 
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of attack of 80 and 50 degrees. Due to the simplicity and speed (short compu- 
tation time) afforded through the use of Eq. (D. 12), this equation was used in 
this study to obtain heating rates off the stagnation lines. 

D. 3 BOUNDARY LAYER TRANSITION CRITERIA 

The following sketches denote the plume impingement phenomena 
occurring during the staging maneuver. 
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Noted in the sketches are several possible shock intersection regions. 

Besides the nozzle lip shock intersection with the bow shock there is an 
additional shock intersection region noted in sketch (B), i.e., the fuselage 
shock intersection with the vertical tail bow shock. 

Several experimental investigations of the leading-edge shock-impingement 
problem have been conducted over the years (Refs. D. 1 1 through D.22). The 
available data for boundary layer transition on the leading edge in the pre- 
sence of an impinging shock wave indicate that transition occurs at Reynolds 
numbers, based on freestream conditions and leading edge diameter, between 
130,000 and 200,000 for sweepback angles of 25 degrees and greater. If the 
freestream turbulence level present in the plume for this study was also 
accounted for, the transition Reynolds number would be lower (Refs. D.23 and 
D. 24. However, in the absence of a detailed study or experimental data on 
freestream turbulence in the orbiter exhaust plume, the lower limit of 
130,000 was used in this study as the boundary layer transition criteria. 

D.4 STRUCTURAL THERMAL RESPONSE 

The structural temperatures presented in this study were calculated 
using one -dimensional transient heat conduction theory as programmed in 
Ref. D.25. Briefly, the material is divided into a series of nodes and cn 
energy balance is performed across each node. For an internal node the 
energy balance equation is 

✓ 



dT 

dt 


k 


dT 

dx 


(D. 13) 


The term or the left hand side accounts for the change in node tempera- 
ture with time and the term on the right-hand side accounts for the transfer 
of heat by conduction. 


s 



ifp- 
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The energy balance for the surface node contains the additional terms of 
heat transfer by convection and radiation and is given by, 


pc 


P 


dT 

dt 


‘conv 


- q 


rad 


+ k 


dT 

dx 


(D. 14) 


The above equations were solved using forward-finite difference tech- 
niques. Also accounted for were the temperature dependency of the material 
prope rties. 
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Fig. D-2 - Chordwise Laminar Heating Distribution 
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